Inflation without Inflaton(s) 



Scott WatsoiE 

Physics Department, University of Toronto, Toronto M5S 1A7 Ontario, Canada. 

Malcolm J. Perry 
Department of Applied Mathematics and Theoretical Physics, 
Centre for Mathematical Sciences, University of Cambridge, 
Wilberforce Road, Cambridge CB3 OWA, England. 

Gordon L. Kane and Fred C. Adams 
Michigan Center for Theoretical Physics University of Michigan, Ann Arbor, Michigan 48109, USA. 

(Dated: May 8, 2008) 

We propose a model for early universe cosmology without the need for fundamental scalar fields. 
Cosmic acceleration and phenomenologically viable reheating of the universe results from a series 
of energy transitions, where during each transition vacuum energy is converted to thermal radia- 
tion. We show that this 'cascading universe' can lead to successful generation of adiabatic density 
fluctuations and an observable gravity wave spectrum in some cases, where in the simplest case it 
reproduces a spectrum similar to slow-roll models of inflation. We also find the model provides a 
reasonable reheating temperature after inflation ends. This type of model may also be relevant for 
addressing the smallness of the vacuum energy today. 

PACS numbers: 



I. INTRODUCTION 

Despite the simplicity and promising phe- 
nomenology of scalar driven, slow-roll inflation, 
much remains to make the idea theoretically vi- 
able. In particular, vexing issues such as the re- 
quired flatness of the inflationary potential and the 
very existence of a fundamental scalar (which must 
be both extremely light and weakly interacting) 
remain elusive (see however [l|, 0, |jj). In recent 
years, a substantial effort has been invested in un- 
derstanding how to embed such models in a quan- 
tum theory of gravity @, 0], and there has also 
been the suggestion of removing the need for slow- 
roll completely Q (see @ for earlier work). How- 
ever, in this paper we will take a different and yet 
complimentary approach to inflation model build- 
ing based on fundamental scalars. 

As the universe expands and cools the fields 
and particles affecting the expansion pass through 
a number of different phases. In the very early 
universe many of these transitions may have been 
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inflationary. In this paper, we propose the idea 
of a 'cascading universe' by asking: if the uni- 
verse passed through enough of these transitions, 
could this provide an adequate alternative to scalar 
driven inflation? Our goal will be to examine 
whether the proposed cascading universe model 
can satisfy the rigid constraints required of success- 
ful inflation model building. We will postpone the 
very important question of embedding the model 
in a fundamental theory (such as string theory) for 
future work. 

In Section [HI we present the cascading model 
and obtain the constraints on the decay rate in 
order for adequate inflation to solve the standard 
cosmological problems (e.g. horizon and flatness 
problems). In Section lllli we address the nature 
of the transitions and we find that adequate in- 
flation and successful reheating are possible given 
certain constraints on the nature of the transitions. 
In particular, we find that cascading can proceed 
via second order or weakly first order phase tran- 
sitions and at a rate r < H as might have been 
anticipated from intuition coming from the grace- 
ful exit problem in old inflation. Working under 
these assumptions, in Section [D71 we turn to the is- 
sue of cosmological perturbations and demonstrate 
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that a nearly scale invariant spectrum of both den- 
sity and tensor perturbations results in the sim- 
plest case of a constant decay rate. We also find 
that in the more realistic case of a varying decay 
rate it may be possible to distinguish this model 
from the usual slow-roll models in that the evolving 
adiabatic sound speed can result in an observable 
tensor to scalar ratio. We conclude in Section [V] 
where we summarize our results and discuss future 
considerations. 



II. A CASCADING UNIVERSE 

Let us consider the case of a universe domi- 
nated by vacuum energy and an additional radia- 
tion component that is subdominant. For a homo- 
geneous and isotropic universe the Einstein equa- 
tions are 
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where we work with the Planck mass, which is re- 
lated to the Newton constant by Gn — M~ 2 . The 
continuity equation is given by 



V M T^ = 0, 



p = -3H(p + p), 



(2) 
(3) 



We consider a two component fluid composed of 
radiation and vacuum energy. The energy density 
and pressure are given by 



P = PA + Pr, 
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with A > 0. For these sources the equations of 



motion (TT|) and ^ become 
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where we have used p\ = 0. We see that in or- 
der for acceleration to occur we need p r < p\- In 
fact, the amount of radiation present is a measure 
of how far the universe is from an exactly de Sitter 
phase. Quantitatively this can be seen by consid- 
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where e is a parameter measuring the deviation 
from a pure dS space-time 1 . For inflation to occur 
we thus expect e <C 1, which for this background 
gives 



2p r 
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One can solve the background equations ([7]) in 
the absence of a coupling and we find 
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are constants chosen so that when t = we have 
a = 1. 



1 This is analogous to the slow-roll parameter e in models 
of scalar field inflation, but because our model does not 
contain any scalar fields we will avoid this terminology. 
Moreover, in contrast to the slow-roll case, the definition 
of e is exact and does not depend on any approximation. 
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Since this solution is derived for the case of A = 
constant, inflation of course does not end, and this 
solution does not provide for a successful inflation- 
ary epoch. When we include time varying values 
of the vacuum energy, which results from multi- 
ple transitions or cascades, we will see the result 
becomes satisfactory. 

The energy transfer from the vacuum energy 
density to radiation (massless string states) is 
given by 



Qa = -I>A, 
Qr = I>A- 



(13) 



where T is the transition rate and the modified 
continuity equation becomes 

V»Tf - pa-Qa, (14) 
V^Xf 5 = p r + AHp r = Q r , (15) 
T^)=0. (16) 
(17) 



The equations of motion for the background are 
then 



3H2 = W p {PA + Pr) ' 



(18) 



PA = -I>A, (19) 

p r = -AHp r + Tp A . (20) 
From (|19|) we immediately find 

PA = PA e- rt . (21) 
Using this result in the above equations we find 
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The solutions are related to modified Bessel func- 
tions. They can be simply expressed by introduc- 
ing the dimensionless quantity 
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The scale factor and Hubble parameter are then 
given by 



= ^ {otiIa{r) + a 2 K (T)) , 



(24) 



= rr / a 2 K 1 (T) - aihjjy . 
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FIG. 1: The deformation parameter for various values 
of the decay rate V. Time is measured in units of the 
Planck time and we take Mp S> pA 2> pr- The various 
curves are given by the values F — 0.01, 0.02, 0.03, 0.05 
from left to right. As discussed in the text, e is initially 
small and proportional to the density of radiation p r . 
As inflation proceeds, more and more energy is dumped 
into radiation via the coupling V. At the very end of in- 
flation we are left with a radiation dominated universe 
corresponding to e = 2. 



where the functions I v and K v are modified Bessel 
functions of order v (see e.g. [35|]). The constants 
are given by 



rv, 

4 

rr, 



-Jfi(T e ) - H e K (r e ), (26) 



a 2 = -^Ii(r e )+H e I ( Te ), (27) 

with r e = r(t e ), H e the Hubble parameter at the 
end of inflation (t — t e ) and we normalize so that 
the number of efoldings is measured from the end 
of inflation, N = ln(a e /a) = — ln(a) where a e = 1. 

We can again introduce a deformation parame- 
ter as in ([5]), however now it is time dependent, 



i(t) = 2pr = 2n r . 

PA + Pr 

Using this in (fT5|) we find 

e(t) = 2- W7TPAo e- rt 



(28) 



(29) 



with H being given by (|25[) . At the beginning of 
I 2 



inflation we have H 2 M 2 ~ p A and so e < 1. As 
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FIG. 2: The number of e-folds of inflation N = / i/dt 
for various values of the decay rate T. Time is mea- 
sured in units of the Planck time and we take Mp ^> 
PA S> Pr- The various curves are given by the val- 
ues T = 0.01,0.02,0.03,0.05 from top to bottom. We 
see that the requirement of sufficient inflation places a 
constraint Y < 0.02. 



the energy is transferred from the vacuum energy 
density to radiation via particle production, the 
deformation parameter increases as can be seen in 
Fig. [TJ Inflation then ends when p r ps p\ and 
e ps 1, which can be seen in Figures [1] and [3] The 
final result is e = 2 and we are left with a universe 
filled by the radiation 2 p r . Let us now consider the 
amount of inflation or number of e-folds. Using 
the above expression for the deformation parame- 
ter the Hubble equation can be rewritten as 



3H 2 = 16ttpa r _ rt/2 



M|(2-e) 



(30) 



This can then be integrated to find the number of 



2 Much later, of course, the radiation is diluted as the vol- 
ume increases and the small remaining constant energy 
density again dominates. 





FIG. 3: The graphs above show the evolution of the 
vacuum energy density SI a = Px/p and the radiation 
energy density Q r = p r /p relative to the total density. 
We present the evolution for two values of the coupling 
r = 0.05 (top) and T — 0.01 (bottom), where it an be 
seen that stronger coupling means inflation ends faster, 
through faster dissipation. 
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where to = is the beginning of inflation, t e is 
the end and is the initial Hubble scale. In the 
second line we use the fact that the denominator 
varies smoothly from y2 to 1 and exp(— Tt e /2) w 
0. As we may have anticipated the amount of in- 
flation depends on the initial vacuum density and 
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the decay rate. 

As an example, consider inflation with a Hubble 
scale near the GUT scale H b ~ M GUT « 10 15 GeV. 
We see to get adequate inflation the decay rate 
need only be slightly below the initial Hubble scale 
r oj 10 GeV . This condition is required in order 
that cascading lasts long enough so that the cosmic 
acceleration can resolve the horizon and flatness 
problems. In Figures 1-3, we examine the evolution 
numerically and find adequate inflation is possible 
given modest values of the parameters. 

Another important consideration is the reheat 
temperature of the model. The cosmic accelera- 
tion ends at the moment t r when p r = p\ and 
radiation comes to dominate. At this moment we 
have 3H 2 = 16irpA where p\ = p r = p^ Q e~ Ttr . 
Using the exact solution (j24|) and (j25|) and assum- 
ing the minimal amount of efoldings (N = 60) we 
find Tt r « 10 so that the reheating temperature 
can be approximated as 

T ~ „i/4 _ 1/4 -rt r /4 

sa 10 15 GeU, (32) 

where we have used pa — An-Mp/87r and we have 
taken the initial Hubble scale Hi, w 10 14 GeV. 
We will see in the next section that this is con- 
sistent with producing the observed temperature 
anisotropics in the cosmic microwave background 
and avoiding over production of gravity waves. 

III. NATURE OF THE TRANSITIONS 

Another important constraint on the cascading 
model comes from considering the type of phase 
transition from level to level. In arriving at the 
previous constraint on the decay rate T in the 
last section, we have tacitly assumed that what- 
ever the nature of the phase transition that it was 
successfully completed and that the associated mi- 
crophysics was irrelevant. However, obviously this 
is not always the case and properties of the transi- 
tions, such as whether they are first or second order 
can play an important role. We will now consider 
both cases of first and second order transitions. 

First order transitions typically proceed by nu- 
cleation of bubbles of the new phase in the back- 
ground of the old phase (see however [27[). The en- 
ergy difference between the phases is stored in the 



bubble walls, and typical expand near the speed of 
light. The transition from the old phase to the new 
phase is complete when all the nucleated bubbles 
collide, releasing the energy stored in their walls. 

Completion of such transitions in cosmological 
backgrounds can often be problematic. If the tun- 
neling barrier (more precisely, the tunneling ac- 
tion) is large, then this corresponds to a strongly 
first order transition. In such cases the average 
bubble size is typically comparable to the gravi- 
tational scale (-ff -1 ) and the gravitational back- 
ground can have important effects [28|, |29|, |30( . In 
particular, in the case of inflationary backgrounds, 
bubbles of the new phase will form in the expo- 
nentially expanding background of the old phase. 
In this case, although the nucleated bubbles ex- 
pand at the speed of light, the background itself 
is expanding faster. This makes bubble collisions 
rare and instead of completion of the phase transi- 
tion we find isolated bubbles expanding in a back- 
ground of eternal inflation. This is the graceful 
exit problem. 

This was made more precise by Guth and Wein- 
berg in [28j | . Consider the zero temperature bubble 
nucleation rate per volume 3 

t = Ae- SE , (33) 

where Se is the Euclidean action and A comes 
from a one-loop determinant factor that depends 
on the microphysics and is typically the energy 
(density) scale of the transition. We note with 
hindsight that if the transition is going to com- 
plete then the typical bubble size r& must be much 
less than the gravitational scale 4 , i.e. ri, -C if -1 . 
In this case gravitational effects in are neg- 
ligible. The authors of 28] then showed that the 
probability of a point to remain in the false vacuum 
is given by 

p(t)~exp(-^/3Ht), (34) 



3 Note this is not the same as the decay rate T that we 
have introduced above, however we will see that these 
quantities can be related. 

4 A detailed discussion of this point appears in [3d| , where 
it is shown that the only feasible first order transition is 
one that results in a distribution of bubble sizes sharply 
peaked around r;, and with -C H~ 1 , i.e. a weakly first 
order transition. 
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where we introduce the dimensionless quantity [3 = 
T/H 4 . We see that the corresponding decay time 
is given by 

The number of efoldings resulting from the tran- 
sition is N w Ht. The authors of [H, H3] found 
that for percolation to occur and the transition to 
complete f3 > (3 C — 9 /(An), which we see corre- 
sponds to N < 1/3. Thus, we see in the cascad- 
ing model we have presented in this paper, that 
if we assume transitions are first order than they 
must be weakly so and the maximum number of 
efoldings per level is N — 1/3 corresponding to 
t/H 4 = 9/(An). To build a successful inflation 
model we will need a large number of transitions 
( w 150 — 180) of this type to get the required 
50 — 60 efoldings in order to solve the horizon and 
flatness problems. This idea for the case of a fun- 
damental scalar field has been argued for in 0, [j| 
and is similar to the original proposal of Abbott 
[3lj |. To relate this to the phenomenological cas- 
cading model we presented in the last section, we 
note that this is exactly what one expects from a 
course grained approach. The decay of vacuum en- 
ergy to radiation will proceed as explained above, 
but on gravitational scales t > H^ 1 this will sim- 
ply be seen as a decay of the vacuum energy to 
radiation, whereas on very small scales t <C H^ 1 
one would see that it was small bubble collisions 
that had been responsible for replenishing the ra- 
diation bath. Thus, by treating the percolation 
events as instantaneous (as seen from the gravi- 
tation scale) the constraint on T (as opposed to 
r) is simply that of the last section T/H < 1/8. 
Combining both the microscopic constraints, along 
with the macroscopic requirements we see that the 
first order transition can work, but it is in a very 
special regime. It will be a challenge to embed 
such an approach into a fundamental theory, such 
as string theory, where e.g. in the case of scalars, 
potential barriers typically lead to strongly first 
order transitions. 

Of course the other possibility is that the transi- 
tions proceed via a second order phase transition. 
Since such transitions do not proceed via quan- 
tum mechanical tunneling, bubble percolation is 
no longer a concern. Immediate examples of such 



transitions are provided by that of a scalar field, 
in which case we simply recover slow-roll inflation 
models such as new inflation and chaotic inflation. 
Another possible example is provided by that of 
the 'thermalons' [32j |. where transitions are stimu- 
lated by thermal fluctuations and so-called 'over- 
barrier' tunneling. 

In summary, we find that the cosmological phase 
transitions of the cascading model must either be 
second order or weakly first order. 



IV. COSMOLOGICAL PERTURBATIONS 

Let us consider density and tensor fluctuations 
about the background solution ([7|) . We will be pri- 
marily interested in modes that leave the Hubble 
radius 50 — 60 e-folds before the end of inflation, 
since these are the modes that are responsible for 
the CMBR anisotropics observed today. There- 
fore, instead of working with the exact solution ([7]) 
for the study of perturbations, it will often be sim- 
pler to work in the conformal time r\ = — oo . . . 
where dr\ = a~ l dt and with approximate solution 

a(r/) = (- V )- {1+S \ (36) 
n( V ) = — , (37) 

where TL is the conformal Hubble parameter H = 
aH and is related to the deformation parameter by 
e = 1 — |t2- The approximate solution treats the 
deformation parameter as a constant, since its rate 
of change is small (e ~ e), but of course this solu- 
tion must break down towards the end of inflation 
when e ~ 1 . 

We now consider linearized perturbations about 
the background ([7]) and in what follows we will 
adopt the conventions of [Hj]. 

We work in longitudinal gauge with the per- 
turbed line clement 

ds 2 = -(1 + 2<5>)dt 2 + (1 + 2* + h i:i )dx l dx j , (38) 

where the tensor perturbation is traceless and 
transverse (i.e., h\ = V — 0) and can be bro- 
ken into its two polarizations h± . Our background 
contains no anisotropic stress so one finds from the 
Einstein equations $ = \f r . Thus, we have only 
one scalar metric degree of freedom associated with 
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the density perturbation and two tensor metric de- 
grees of freedom for the gravity waves. Because we 
are considering linearized perturbations the scalar 
and tensor metric fluctuations decouple and we will 
treat each one in-turn. 



A. Density Fluctuations 



Working in conformal time the equations for 



sistent with the fact that a true cosmological con- 
stant does not propagate, i.e. is constant. This 
allows us to see the importance of the graviton 
and other particle production (T term), since in 
an quasi-exponentially expanding background with 
no transfer the perturbations will be immediately 
damped away. 

The equation of state parameter is given by 



density perturbations are, 

V 2 $ - m{H<5> + $') 
di (Hi + *') 
i" + 3W$' + (2H' +H 2 )<f> 



4ttGci 2 5p (39) 
4irGa(p + p)5u h (40) 
4irGa 2 6p, (41) 



t = 
P 

= -1 



-1 



(45) 
(46) 



where Sp, 8p, and Sui are the perturbations of the 
total energy density, pressure, and velocity, respec- 



tively and H 



da 
dn 



is the conformal Hubble 



parameter and V is the comoving gradient. The 
pressure perturbation is related to the energy den- 
sity and entropy density perturbations by 



Sp = — — 



c , dp 
os 



dp 
dp 

c 2 Sp + tSS, 



SS 7 



(42) 
(43) 



with c s the adiabatic sound speed at which the 
perturbations evolve and should not be confused 
with the equation of state parameter w = p/p, 
which depends on the background quantities. 

Indeed, for the model we consider here these 
quantities are quite different. The adiabatic sound 
speed for long-wavelength perturbations is given 
by 5 



P 
P 



3Hp r 



(44) 



where we see in the limit r — > our perturba- 
tions will evolve like pure radiation. This is con- 



5 Strictly speaking it is incorrect to think of = ? as the 
sound speed during inflation, since this is only true on 
large scales where (? s = | ~ 4 and these modes evolve 

on scales beyond the sound horizon. On small scales dur- 
ing inflation the metric perturbation $ oscillates and the 
effective adiabatic sound speed is found to be (? s = 1, in 
agreement with causality. 



which reduces to the pure de Sitter solution if 
p r = e = as we have noted. We see that w(t) is 
explicitly time dependent and in general this im- 
plies that there can be a significant contribution 
from non-adiabatic pressure in (|43[) by the entropy 
term tSS. This could result in significant genera- 
tion of entropy perturbations, a possibility that we 
will analyze in Section pV A 2p . 

We now return to solving the system (f3"9"| - (|41[) . 
Combining equation (|39|) with (|4ip and working in 
momentum space V 2 <& — > — /c 2 $, we find a second 
order differential equation, 



+ 3W(l + c 2 )$' fc + 
+ (l + 3c 2 )H 2 ]$, 



2, 2 



ZH'+ (47) 
4irGa 2 TSS. (48) 



subject to the constraint (|4U)) . We can simplify 
this equation by introducing the field redefinition, 



ik 



Uk 



AnGp 1 / 2 ^ 
8ttM~ 2 



w 



3a 2 (l + w)' 
where w =p/p. Then (|47p becomes 



; 2 2 

k c% 



(49) 
(50) 

(51) 



with M giving the contribution from entropy 
modes as 



M = cPp^VT+w tSS. 



Adiabatic Fluctuations 



(52) 



We will first consider solutions to (|51|) in the ab- 
sence of entropy modes, i.e. JV = 0. For modes 
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that are far inside the sound horizon kc s S> Ti we 
find that u oscillates with a constant amplitude 
that is to be found by the initial conditions after 
quantization. Far outside the horizon kc s <C TL 
and by inspection we have the solution u ~ 6. 
However, this solution corresponds to a decaying 
mode for the metric perturbation Instead, it is 
the growing mode that is of interest, which can be 
found by noting 9" w during inflation, so that u 
constant is also a solution. Assuming that we are 
deep in the inflation epoch where e is nearly con- 
stant an exact solution can be found by integra- 
tion. In terms of the original metric perturbation 
$ one finds 



$ = 



U i , 



H f . 2 A 
= aA Ax+ 3 



.!„. / (l + w)a 2 dr/ 
ea 2 drj 



H 



(53) 
(54) 
(55) 



where we have used 



w — — 1 + |e is nearly con- 
stant during inflation. The first term in (|55jl corre- 
sponds to the decaying mode found above (u ~ 9), 
whereas the second mode is nearly constant with 
$o = to be determined by matching to the os- 
cillating mode inside the sound horizon. We note 
the importance of the graviton production in this 
model resulting in a non-zero radiation density 
since in the pure de-Sitter case where p r = so 
that e = we see no density metric perturbation 
would remain since <E> — * as e — * and all that 
is left is the decaying mode. This is an illustration 
of the no hair theorem for pure de Sitter space. 

In summary, we have found that during infla- 
tion the metric perturbation is nearly constant on 
super-horizon scales, whereas on sub-horizon scales 
we find u ~ $ undergoes constant amplitude oscil- 
lations. What remains is to quantize the perturba- 
tions in order to determine the unknown constant 
$o- However, we must first justify neglecting the 
entropy mode term ( i.e. M) in (IB"!]) . 



2. Entropy Flucutations 

In this section we consider the role of entropy 
fluctuations in the model. We will follow [34[ where 



a systematic procedure for the study of perturba- 
tions in multi-fluid systems was described. It will 
be useful to introduce (, which is curvature pertur- 
bation on constant energy density hypersurfaces. 
We will drop the momentum index in what follows, 
writing ( = ( k . 

In the presence of multiple fluids, the total cur- 
vature perturbation can be expressed as a sum of 
the curvature perturbation due to each fluid com- 
ponent as 



where 



1~L 

(a = <f> + —5p a , 



(56) 



(57) 



and we have used the lack of anisotropic stress to 
again write $ = as before. For entropy fluctua- 
tions we are interested in the non-adiabatic contri- 
bution to the pressure perturbation in (|43[) , which 
is given by 



Sp r , 



tSS = Sp — c^Sp. 



(58) 



As discussed in [34| , there are two sources of non- 
adiabatic pressure 



Sp„ 



Sp, 



rel 
nad 



Sp. 



int 
nad 



(59) 



which are the relative and intrinsic non-adiabatic 
pressures, respectively. In the model we are con- 
sidering here the two fluid components have fixed 
equation of state, i.e. Spa = —Spa and Sp r = 
1/S5p r so that there is no intrinsic non-adiabatic 
pressure, i.e. dp™* d = 0. The contribution to the 
relative non-adiabatic pressure is [34| 



a.f3 



i 



-pApr(c\ - cl)SAr, 



(60) 



where we have introduced the relative entropy per- 
turbation 



Sa/3 = 3(Cq - 0), 

= -m ( s Jsl 



<>££ 



(61) 
(62) 



9 



with the factor of three due to the convention of 
normalizing to baryons. For the model we consider 
here 

S Ar = -S rA = -3h( S 4±- 6 -^). (63) 
V Pa Pr J 

Returning to ([60)) we see that the relative non- 
adiabatic pressure is proportional to p r , the rate 
of change of the radiation density. During infla- 
tion one finds from the background solution (|24p 
and ([25)) that p r w 0. That is, the transfer of vac- 
uum energy to the radiation density via the cou- 
pling r is just enough to counter the dilution of the 
radiation by the exponential expansion. Thus, the 
non-adiabatic pressure is negligible and we need 
not worry about the presence of entropy perturba- 
tions during inflation. 

However, there is a more fundamental reason to 
expect entropy perturbations to be absent from 
this model. The crucial point is that a relative 
entropy perturbation is produced when two fluids 
generate different curvature perturbations. This 
difference can then be mediated from one fluid to 
the other via the gravitational background. A well 
known example is the perturbation in the baryon- 
photon ratio 

Sb 7 = 3(Cb - C 7 = T — L , 64 

PB 4 p 7 

which does not vanish because the two fluids are 
perturbed differently. 

However, in the case we consider here things are 
different. In the absence of the coupling F there 
is only one fluid with propagating fluctuations, 
namely the radiation density with fluctuations Sp r . 
In this case the long-wave fluctuations propagate 
at Cg = 1 /3 and the cosmological constant remains 
a constant, i.e. 5p\ — 0. In the presence of the 
coupling r the fluctuations now propagate at a dif- 
ferent adiabatic sound speed , but the two flu- 
ids are coupled through their equations (|19p and 
([20)) through the term ±Tp\. Thus, there is really 
only one propagating degree of freedom and the 
two fluids do not evolve independently, resulting 
in S\r = 0. 

In fact, in this regard this is not unlike the 
case of inflation by a single scalar field where it 
is known that there are only adiabatic pertur- 
bations. Instead of working with the scalar di- 



rectly, we could consider two fluids, one represent- 
ing the kinetic energy with a stiff equation of state 
pi = pi = l/2cb 2 and a second fluid composed of 
the potential pi — —pi — V{4>). Insisting on this 
two-fluid description and demanding that the full 
equations of motion are satisfied we are led to an 
energy exchange term Q± = ±(j)V'(<f>), similar to 
the case we have above. However, since we know 
there is only one degree of freedom, we certainly 
know that there are no entropy perturbations and 
no non-adiabatic pressures. This can be seen by 
examining the perturbation equations in full detail, 
and in particular one finds that the two fluids do 
not evolve independently due to the coupling Q± 
and the fact that the second fluid does not prop- 
agate in the absence of the coupling (i.e. 5p2 = 
for Q = 0). 

In sum, we see that entropy perturbations in the 
cascading model are negligible during inflation for 
the case of a constant decay rate F. For the case 
of a time varying T, this issue must be revisited, 
which is work in progress. 

B. Spectrum of Fluctuations 

Having shown that entropy perturbations are 
negligible, we proceed to find the spectrum of the 
density fluctuations. In order to find the power 
spectrum all that remains is to determine the un- 
known constant <f>o in ([53)1 . We can then find the 
gauge invariant, comoving curvature perturbation 

which is related to the curvature perturbation 
from the last section by 

1 fc 2 $ fc 

** = &+3W = 7F' ( 66 ) 

so that for large scales modes (which are the ones 
of interest) k — > and TZk — ► Cfe- The density 
power spectrum is then defined as 

n = TTlCfc| a , (67) 

which can be compared with observations. 

Finding the constant $ is accomplished by en- 
forcing the correct initial condition on the modes. 
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However, these modes are born in their vacuum 
state far below the Hubble radius. This requires us 
to quantize the perturbations, starting their evo- 
lution in the standard adiabatic vacuum. Then we 
have seen that the solution inside (kc s ^S> Ti.) oscil- 
lates with constant amplitude until Hubble radius 
crossing where it can be matched to the solution 
outside (fcc s -C Ti) providing us with the required 
normalization constant. 

The only obstacle to quantization is finding the 
canonical field which diagonalizes the action. For 
the case of hydrodynamical fluids, as we consider 
here, this was done in [33]. There it was found 
that the canonical field Vk (the so-called Mukhanov 
variable) which is related to Uk by 



Uk 



(v k 0)' 



(68) 



and the curvature perturbation (|65[) by Vk = 
reduces the action to that of a harmonic oscillator 
with time dependent frequency. In terms of this 
variable the equation of motion (|51[) becomes 



z 
z 



v k = 0, 



(69) 



where z = (c s &)~ x . In terms of the new variable Vj~ 
the solutions on large scales (kc s ^S> Ti) are given 
by Vk ~ z. Notice this is the growing mode of inter- 
est in contrast to the classical case where Uk ~ 9 
decayed and it is this squeezing of the quantum 
state that will result in classical fluctuations on 
large scales. On small scales the momentum term 
dominates and we again have oscillations with con- 
stant amplitude. 

We could now proceed with the approximate so- 
lution, however in the case e <C 1 we can solve (j69]) 

exactly. We find — = - ~l' 4 where v — 3/2 + e 
and the solutions can be expressed in terms of Han- 
kel functions. We require that the modes begin in 
the adiabatic vacuum, which amounts to the con- 
dition 



v k 



1 



-ikc s 7j 



as kc s rj — > — ao. (70) 



\/2c s k 

The appropriate solution is then given by 

ot / 1 \ 

v k{V) = 1pJ-f]n\, >(-kc s Tl), 



(71) 



where is a Hankel function of the first kind 
and | a | = 1. We can then immediately find the 
curvature perturbation 



l^jHP(-kc sV ) (72) 



On large scales using the asymptotic expansion of 
the Hankel function we have 



la 



i 



{-kc s fj) 1/2 v , 



(73) 



z^Jlixkcs 

and using z = (c s #) _1 the power spectrum (|67p is 



Pt 



1 



H 



An 2 c s e IM 



(-kc s rj)~ 



(74) 



We note that this reduces to the standard slow-roll 
inflation result for the case \c s \ = 1. 

The tilt of the power spectrum is given by 



n, = 1 



dln"P c 
G?ln fc 
= 1 - 4fi r , 



= 1 - 2e, 



(75) 
(76) 



where f2 r = p r j p. By noting that p r w during 
the time modes of interest exit the Hubble radius 
(i.e. N ~ 50), we see that the tilt of the spectrum 
is set by the initial abundance of radiation since 
e sa 2£l ro is constant during inflation. 

Comparing ^ to the best fit WMAP3 data 



Pr = 19.9 



+ 1.3 



x 10 



-10 



we find that 



H 



1/2 



< 10 _5 M K 



(77) 



(78) 



Since e«l during inflation this implies an upper 
bound on the Hubble scale during inflation H < 
10 14 GeV. Combining this with the constraint for 
adequate inflation from (f3"Tj) . i.e. T/H < 1/N we 
find an upper bound on the decay rate of the vac- 
uum energy T < 10 13 GeV consistent with our 
earlier results and our general approach. 



1. Gravity Waves 

The gravitational wave spectrum is found in 
much the same way as the spectrum of density 
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perturbations. One first decomposes the gravi- 



ton into its two polarizations and hY '. The 
modes then obey the same equation (|69|) as the 
density fluctuations except in this case we have 
Vk = ahkMp (where we set h = h±) and z is re- 
placed by the scale factor a. The solution is again 
in terms of Hankel functions, and one finds for the 
long-wavelegnth fluctuations 



,(-) 



(-kn) 



1/2-1/ 



and the power spectrum is 



Ph = 



H 



{-kr,) 



-2c 



(79) 



(80) 



with the tilt of the tensor spectrum Ut = — 2e. 
Thus, we see that the main difference between 
the tensor and density spectrum is the deforma- 
tion factor and the presence of the adiabatic sound 
speed in the spectrum of density perturbations. 
Our tensor to 'scalar' ratio is then 



Ph 1R . 
— — = loec s 
Pr 



Scalar Free Model (81) 



which contains the adiabatic sound speed c s eval- 
uated at the time of Hubble radius crossing. 

This is an important result and is similar to 
models of kinetic inflation Q , where the adiabatic 
sound speed offers a way to distinguish this model 
from standard slow-roll inflation which gives in- 
stead 



16e, Scalar Slow-roll Inflation, 



(82) 



where we recall that e ss e is the usual slow-roll 
parameter which measures the slope of the scalar 
field potential in units of the Hubble scale. 

We have seen the adiabatic sound speed does not 
differ greatly from the usual slow-roll inflation case 
for the choice of T = constant that we have consid- 
ered here. However, for the case of non-constant 
r this could dramatically change, since the adia- 
batic sound speed could differ greatly from one. 
This could allow for an observable tensor to scalar 
ratio, where standard models of scalar driven in- 
flation starting near the string scale seem to gener- 
ically predict an unobservable spectrum [38[ . This 
is work in progress. 



log H 

(GeV) 



Allowed Region 



N>60 



/ -22 -20 
T > (Exit Inflation) 



5 10 14 l og r (GeV) 



FIG. 4: The figure above summarizes the constraints 
found on the Hubble scale during inflation and transi- 
tion rate V between levels. In addition to these con- 
straints one should also add the requirement that the 
transitions proceed via second order or weakly first or- 
der phase transitions. 



FURTHER CONSIDERATIONS AND 
CONCLUSIONS 



In this paper we have considered a cascading 
model for the early universe that provides a period 
of cosmological acceleration, which can account for 
the required number of efoldings. As the universe 
cascades, vacuum energy is converted into radia- 
tion inhomogeneously, resulting in a nearly scale 
invariant spectrum of cosmological density pertur- 
bations and a small amount of gravitational waves. 
Once the radiation density overtakes the decay- 
ing vacuum energy, the model naturally exits in a 
radiation dominated universe with a temperature 
which we found can be as large as T r w 10 15 GeV. 
As the universe evolves through the radiation and 
matter epochs the vacuum density will once again 
dominate the energy density if the decay does not 
proceed to zero vacuum energy. 

We have seen that our approach has one basic 
(in principle calculable) parameter, the level decay 
rate T. The number of e-foldings, the reheating, 
and the density fluctuations all depend on T, and 
we find there does exist a range of values of T con- 
sistent with the data for all of these, which might 
not have happened. 
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Although these preliminary findings are promis- 
ing, much remains to be addressed. A particularly 
pressing issue is a concrete derivation of the decay 
rate T or equivalently a better understanding of the 
level spacing and the time spent in a given energy 
(density) level. In fact, we argued in Section |V] 
that if r is not taken constant, the result is a vary- 
ing adiabatic sound speed which can result in den- 
sity perturbations and gravity waves that would 
further distinguish the cascading model presented 
here from usual slow-roll inflation. 
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